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STUDY OF THE OSCILLATORY MOTION OF MANNED VEHICLES 
ENTERING THE EARTH’S ATMOSPHERE 
By Simon C. Sommer and Murray Tobak 


SUMMARY 


An analysis is made of the oscillatory motion of vehicles which 
traverse arbitrarily prescribed trajectories through the atmosphere. 
Expressions for the oscillatory motion are derived as continuous functions 
of the properties of the trajectory . 

Results are applied to a study of the oscillatory behavior of re-entry 
vehicles which have decelerations that remain within limits of human toler- 
ance. It is found that a deficiency of aerodynamic damping for such 
vehicles may have more serious consequences than it does for comparable 
ballistic missiles. 


INTRODUCTION 


Studies of the oscillatory behavior of missiles entering the atmos- 
phere on ballistic trajectories (refs. 1 and 2) have revealed that the 
rapid increase in atmospheric density experienced by such vehicles is a 
potent factor in restraining the magnitude of their oscillations. In 
reference 2, for example, it is shown that even a vehicle that is dynami- 
cally unstable in the usual aerodynamic sense will undergo oscillations 
that are convergent over either all or the major portion of its path 
through the atmosphere. Generally, the conclusion has been that for 
ballistic missiles, dynamic instability is a potential source of difficulty 
only in their terminal phase of flight . 

For the case of descending manned vehicles, the situation is not yet 
as clear. Since manned vehicles cannot be permitted to develop the very 
large decelerations experienced by ballistic missiles, means of reducing 
the decelerations must be introduced. A study of this problem by Chapman 
(ref. 3) has revealed that decelerations can be held within tolerable 
limits by the use of small amounts of lift and by the use of re-entry 
trajectories starting from very small initial flight-path angles. The 
question arises: What effect do these efforts to reduce deceleration 

have on the oscillatory motion? 
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A study of this question was undertaken by the authors, using as a 
starting point the results of an analysis of oscillatory motions over 
arbitrary trajectories presented by Tobak and Allen in reference 4. In 
the course of this study, it was discovered that a number of limitations 
and assumptions of the previous analysis could be eliminated, thus both 
widening its scope and simplifying its use. The purpose of the present 
paper is therefore twofold: first to re-examine the motion analysis of 
reference 4 to show how it can be extended, and second to apply the 
results to a study of the oscillatory behavior of a class of vehicles 
which have decelerations within the limits of human tolerance . 


SYMBOLS 


reference area 


-D 


Cl 


■ y m 


drag 


drag coefficient, 
lift coefficient, ^ 


qA 

lift 


pitching moment 
pitching-moment coefficient, 


cl, 


a 


■^a 


Sx 


■'mq 


-m 5 


rate of change of lift coefficient with angle of attack, 


da 


a-+o 


rate of change of pitching -moment coefficient with angle of 


attack, 


bC 


a-* o 


rate of change of moment coefficient with time rate of change 


, a l 

of angle -of -attack parameter *77" , 


ac m 

ddl/V 


a-*- o 


rate of change of moment coefficient with pitching velocity 

, Ql f 

parameter 


v ' \dQi/vJ ^ 0 


rate of change of moment coefficient with control deflection. 


d5 


b-+o 


D 


drag 
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g 

I 

Jo( ) 

K 

L 

l 

M(s) 

m 


parameters defined by equation (6) 
acceleration due to gravity 

pitching moment of inertia about center of gravity 
Bessel function of first kind of zero order 
dynamic stability parameter (eq. (9)) 
lift 

body length and reference length for moment coefficient 
evaluation 

parameter defined by equation (12) 
vehicle mass 


Pi(s)l 

P 2 (s)J 

1 


parameters defined by equation ( 9 ) 
dynamic pressure, pV 2 


r distance from center of earth to vehicle 


s 


distance measured along path of "static” trajectory (see text) 


t time 

u horizontal component of flight velocity (sketch (a)) 

u ratio of horizontal component of flight velocity to satellite 

velocity. - 7 = 
vgr 

V flight velocity (sketch (a)) 

- V 

V ratio of flight velocity to satellite velocity, ~j= 

v vertical component of flight velocity (sketch (a)) 

W vehicle weight 

y altitude 

Y 0 ( ) Bessel function of second kind of zero order 

X,Z axes fixed in space with origin at center of earth (sketch (a)) 
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a 


angle of attack (sketch (a)) 


0 


density parameter (eq. ( 32 )) 


7 


flight -path angle (sketch (a)) 


5 


control deflection angle 


e 


angle of pitch measured from axis fixed in space (sketch (a)) 

9 


angle of pitch measured from local horizontal (sketch (a)) 

P 


air density 


Po 


air density at sea level 


a 


radius of gyration 


* 


angular displacement of vehicle from fixed space axis, Q 
(sketch (a)) 

- 0 

O 


- ( ) 
dt v ' 


( )’ 


ds ' 


( )i 


initial value 


( )max 

envelope curve of oscillatory motion 


( )o 


oscillatory variable 


( >8 


nonoscillatory variable 




ANALYSIS 


Our purpose in the ensuing analysis is to show how a further improve- 
ment can he made in the expression developed in reference k for oscilla- 
tory motions over arbitrary trajectories. Specifically, we shall elimi- 
nate both the assumption of constant aerodynamic coefficients and the 
necessity of breaking the trajectory into straight-line segments. 



Equations of Motion 


The equations of motion defining the vehicle’s path and its 
oscillations about that path may be written as 
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The angles a, 7 , 0, <f>, 0 are defined in sketch (a). Note in both 
sketch (a) and equations (l) that we have taken this opportunity to cor- 
rect an error that appears in the derivations of references 1 and 4; 
namely, the erroneous use therein of a moving rather than.a fixed axis 
to define angular velocity ( 0 ) and angular acceleration ( 0 ). This error 
does not affect the equations for oscillatory motion given in references 1 
and 4, however, as will be seen. 



Sketch (a) 
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Following the argument of reference ve now separate equations (l) 
into two sets of equations, one set (7 S , a s , 0 S ) defining the "static 
trajectory" of the center of gravity of the vehicle, the other ( y 0 , Oq, 0 o ) 
defining the oscillatory motions of the vehicle about the static trajectory. 
Thus, under the restriction |7 0 /7 S | « 1, we get for the static trajectory 
equations 

I 

-mV - C D qA + mg sin 7 S = 0 


r 


mV 7 s + qACp a s + m^-j — g Jcos 7 S = 0 

(X v 


( 2 ) 


a s Z 


0 S Z 


I©s 1 CiuqOLs + Cn^ y + C m ^ y 


m„ v + C m & & 


0 


and. for the oscillatory motion 

mV7 0 + qAC L a Q = 0 

•• f a 0 Z 0 q1\ 

IBo “ qAZ( C^Oq + y” + y J = 0 


(3) 


We note, however, that = u/r, which is a nonoscillatory quantity. 
Hence , we may write 


®s = e s “ ^ 
®0 ” °o 


oo 


so that equations (3) for the oscillatory motion agree with equations (5) 
of reference 4. Again, use of the equality 0 Q = a Q - 7 0 permits the 
set of equations (3) to he combined into a single equation for the 
oscillatory angle of attack: 


with 


Oo(t) + fi(t)oo(t) + f 2 (t)ao(t) = 0 (5) 


fx(t) 


f 2 (t) 


c £™/ c +c ll 

C L a 2m ymg_ + Sy 21 

pV 2 AZ d ( pVA\ C m q c L a / P VAZ\ 2 

"Sx, 21 + dt ^ C L a 2my " Im V 2 / 






( 6 ) 


Equations (2) through (6) retain with reference 4 the implicit assumption 
that the drag coefficient Cp is independent of angle of attack and 
pitching velocity. However, it need not be assumed that the aerodynamic 
coefficients in equations (2) through (6) are independent of Mach number. 


7 


Oscillatory Motion 


As in reference 4, henceforth the subscripts s and. o will be 
omitted. It is to be understood that in referring to a we mean the 
oscillatory angle of attack, whereas in referring to 7 we mean the 
flight-path angle of the static trajectory. Now consider equation (5) 
for the oscillatory angle of attack and let the independent variable be 
distance traveled along the flight path, s. The use of s rather than 
y as independent variable has a twofold advantage: First, s is a 

single -valued function for any flight path and, second, the appearance of 
sin 7 is suppressed. With the substitutions 


a(t) = a' (s)V 

a(t) = a" ( s ) V 2 + a' (s)W'(s) 


(7) 


in equation ( 5 ), we get for the oscillatory equation of motion 


a" (s ) + a' (s ) 


V (s ) 


- ? x ( b ) 


+ P (s)a(s) = 0 


(8) 


with 



(9) 


and where we have retained only the dominant term in P 2 (s) (cf. appendix 
of ref. 4). The transformation 


then gives 


<x(s) 



rV’(s) 

V(s) 



ds 


a" (s ) + M(s)a(s) = 0 


( 10 ) 


(H) 


with 


M(s) = P 2 (s) 


1 _d r v(s) 

2 ds V 



1 rv>(s) 
4 V 


P x (b) 


-t2 


( 12 ) 
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Again as in reference k, it is easy to show that M(s) is essentially 
P 2 (s). 


Asymptotic solution for a(s) .- In reference k, it was found that 
with the assumption of constant flight-path angle and constant C m a 

solution to equation (ll) could be written in terms of zero-order Bessel 
functions of the first and second kind. Experience with this solution 
has shown^ however; that it approaches its asymptotic form very rapidly. 
This suggests that the asymptotic solution to equation (ll) likewise 
should be of sufficient accuracy for most practical purposes. The advan- 
tage to seeking an asymptotic solution to equation (ll) is of course that 
the above-mentioned assumptions need not be invoked. It is by this means; 
therefore; that the necessity of breaking the trajectory into a series of 
straight-line segments is eliminated. 


To find the asymptotic solution to equation (ll) we use a theorem 
given in its most general form by Wintner in reference 5* Wintner has 
shown that for large s any differential equation of the form (ll) whose 
coefficient M(s) satisfies the conditions 


M(s) > 0 for all s 


(13) 


00 



[M(s ) ] 


1/2 

5 

fM't 

M" 


4 

W 

“ W 


ds < 00 


(HO 


will have a general 
the form 

a(s) 


solution which approaches asymptotically (as 

“ 1/4 

= [M(s)] [c-jcos cp(s) + c 2 sin cp(s)] 


S -* co) 


(15) 


where 


cp(s) 



(16) 


Further; differentiation of the asymptotic representation (15) is allowed. 
Since; as has been mentioned; M(s) is essentially equal to P 2 (s); the 
condition M(s) > 0 will be satisfied so long as the vehicle is statically 
stable (C^ < 0 ) and its flight path remains within the measurable atmos- 
phere (p > 0). With M(s) > 0 and twice differentiable; the condition ( ik ) 
is readily satisfied for any flight path that terminates in a finite 
distance . 


Hence; letting M(s) equal P 2 (s) combining equation (15) with 
( 10 ), we have as an asymptotic solution for the oscillatory motion 


a(s) = 


1 

1/4 


_1 

2 


P 1 (s)ds 


[V 2 (s)P 2 (s)] 


e 


[cjcos cp(s) + c 2 sin cp(s)] ( 17 ) 
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or 


a(s ) 


[-C^CsJqCs )] 174 


X 

2 

e 


Pjsjds 

cos [cp(s ) - cp D ] 


(18) 


with 


cp(s) 



ds 


Because of the mildness of the conditions affecting its generality, 
equation (l 8 ) is considered to be the central result of this analysis* 
However, it cannot be expected to apply accurately over the initial part 
of the motion. On the other hand, the Bessel function solution of refer 
ence k does apply to the initial motion, since over a sufficiently small 
interval sin y and ^m^ will not change significantly# A simple 
artifice that combines these observations to give one expression that is 
applicable over the initial part of the motion and still has the correct 
asymptotic behavior is the following: Replace the asymptotic solution (15) 

by the Bessel function combination that has (15) as its asymptotic 
representation. Thus 


a 


(s) = 


cp 2 ( s ) 

_ M(s)J 


U/4 


c 3 J Q [cp(s)] + c 4 Y 0 [cp(s ) ] 


(19) 


and 


a 


(s) = 


q> 2 (s) 

V 2 (s)M(s)J 


14 2 I P 1 (s)dS 


3Jo[t( S )] + c 4^ot ( ?( s )] 


( 20 ) 


where M(s) ~ P 2 (s) and - 9( s ) continuous variations in Cn^ and y are 

to be permitted. It can be verified that with and sin 7 constant, 

cp2( s ) = cM(s ) , so that equation (20) then reduces to the form given in 

reference 

Further simplifications.- Equation (l 8 ) holds for any flight path 
that remains^ within the measurable atmosphere, and for arbitrary varia- 
tions of the aerodynamic coefficients . When some of the coefficients can 
be considered constants, further simplifications are possible. As an 
example, let us assume that all coefficients are constants. The exponent 
in equation (l 8 ) then has the form 


JL 

4 


CpA N 

m 


K / p(s )ds 


e 


( 21 ) 
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The integral can be evaluated by adopting the approximations introduced 
by Chapman in reference 3 . Thus, consider the equation for balance of 
forces in the horizontal direction (cf. ref. 3): 


m (dt + t) + D cos ^ " d tan = 0 

Following Chapman, we neglect the terms uv/r and (L/D)tan 7 
tion (22) under the conditions that 

| uv/r | | dr/r| 

| du/dt | ~ | du/u| ^ 

^ tan 7 I « 1 


( 22 ) 

in equa- 


( 23 ) 


The first of these conditions has been shown in reference 3 to be satisfied 
by any trajectory that has entered the atmosphere, since as soon as the 
drag becomes important in decelerating the motion the percentage change 
in u becomes large compared with the percentage change in distance from 
the planet center. The second condition is of course satisfied without 
approximation by nonlifting vehicles at any flight-path angle other than 
90 . It should also be easily satisfied by lifting vehicles entering the 
atmosphere from satellite orbits for then tan 7 is generally very small. 
Equation ( 22 ) becomes 


du 1 f CpA \ u 2 

dt 2 m J P cos 7 


= 0 


( 24 ) 


and the transformation du/dt = (u/cos 7)(du/ds) permits p(s) to be 
written as 



1 du 
u ds 


( 25 ) 


Inserting equation ( 25 ) in ( 2 l) and integrating then gives for <x(s) 

-K/2 


( s ) = / ~~~i 74 cos[<p(s) - cp 0 ] 
[q.(s)J 


(26) 


The envelope of the oscillation has the simple form 

-K/2 


“max ( s ) — 


Cju 

tq(s)] 


1/4 


(27) 
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Convergence Criteria 


General criterion .- Let us now return to the asymptotic solution for 
oscillatory motion (eq. (l8)) and attempt to derive from it a criterion 
whose satisfaction ensures that the motion is convergent at any position 
s of the flight path. We use equation (l8) rather than ( 26 ) since it 
turns out that a remarkably simple result can be obtained from equation (l8) 
free of approximations either to the aerodynamic coefficients or to the 
flight path. Thus, noting that the envelope of the oscillation varies as 


cknaxC 5 ) ~ 


Jr J P^sjds 


[ )q_(s ) ] 


1/4 


ve have convergence at any position s provided ^ (amax) ^ 0. 
ensured if 


( 28 ) 


This is 


F i( s) 1 <l'(s) 1 Gmq^s) 

2 ^ q(s) 4 Cmjs) ° ^ 

It will be noted that equation (29) has essentially the same form (for 
Cjjj^ = const.) as the convergence criterion derived in reference k. Equa- 
tion (29) is the more convenient expression, however, in that it applies 
without change to any trajectory. With the use of y rather than s as 
independent variable the inequality of reference !+ must be reversed when- 
ever 7 changes sign. 

The usefulness of the criterion (29) is enhanced if an expression in 
terms of primary quantities can be found for q'(s)/q(s). Such an expres- 
sion is obtainable from the first of the static trajectory equations 
(eq. ( 2 )) in the following way. Converting the variable in equation ( 2 ) 
to altitude y by the transformation 


dY dV d^ 

dt dy dt 


t . w dV 

•(sm 7 )V ^ = 


sin 7 d 
2 dy 


(V 2 ) 


( 30 ) 


ve have 


1 

2 dy 


(V 2 ) = 


•CpA q 

m j sin 7 


g 


(31) 


Now introduce the well-known assumption that atmospheric density varies 
exponential! y with altitude : 

-fry 

p = P Q e 


(32) 
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Use of equation (32) in the identity 


1 iL (y2 ) _ A 

2 dy ' dy 



then gives 


dq /CdA\ q p 

dy _ \ 111 ) sin 7 ^ 


(33) 

(34) 


and 


1 dq 
q ds 


i dq 

-sm 7 — 

q dy 


= -sm 7 


CpA" 
m 


sm 7 


M 


(35) 


Inserting equation (35) and also the definition of P x (s) (eq. (9)) in 
equation ( 29 ), we get finally as a general criterion to be fulfilled by 
the dynamic stability parameter K 


K < - 1 


f— 

)(sin 7 ) 

"l P 

1 c ma'( s ) 

\CdA, 

[_q + gp + gp 

Sin ? Cmjs) j 


(36) 


where 



Approximate criterion .- The derivation leading to equation (36) has 
not required the introduction of approximations either to the trajectory 
equations or to the aerodynamic coefficients. When equation (27) is used 
to represent the oscillation amplitude rather than equation (28), however, 
it is advisable to use a convergence criterion that is consistent with the 
approximations introduced in the derivation of equation (27). With the 
approximations indicated by equations (23) in force and with constant 
aerodynamic coefficients the convergence criterion obtained from differ- 
entiation of equation (27) takes the form 


K u» (s) _ 1 q 1 (s) 
2 u(s) 4 q(s) 


(37) 


Under the conditions defined by equations (23), the derivatives indicated 
in equation (37) ca *i be evaluated from the static trajectory equations 
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used, by Chapman in reference 3 * (The berm u T (s)/u(s) is already avail- 
able as equation (25).) The convergence criterion becomes 


K < - 1 


L sin 27 ( w V • a 

n ~2— + (^A/ Sin 7) 


D 


P 1 

gp + <1 (1 



( 38 ) 


Oddly enough, it turns out that the more precise expression, equation (36), 
is also simpler in form. However, so long as the conditions imposed by 
equations (23) are satisfied, the differences between results computed 
from equations (38) and (36) should be insignificant . 1 

To summarize, it is advised that equations (28), (29), and (36) be 
used when the properties of the trajectory are known precisely. The 
results (27), (37), and (38) are those appropriate for use with constant 
aerodynamic coefficients and in conjunction with trajectories computed 
from Chapman's approximate analysis (ref. 3 )* 


DISCUSSION 


Let us now apply the results just developed to a study of the oscil- 
latory behavior of vehicles whose decelerations during their descent 
through the atmosphere appear to be within the bounds of human tolerance . 
We have chosen as examples a series of nonlifting and lifting vehicles 
with values of entrance angle ranging from 0° to 4 ° and drag-loading 
parameter, W/CpA, of 10 , 30 , and 100 pounds per square foot. For these 
examples the “static flight trajectories were computed by the method of 
reference 3 with the following conditions common to all cases : 

(a) Constant aerodynamic coefficients 

v i 

(b) Vi = = 1 

*/gr 

(c) yp = 400,000 ft 

(a \ " Py 

(d) P = P 0 e 

with 


^-Note that the spurious l/d term in equation (38) discarded 
in comparison with the first term, since by virtue of the second of equa- 
tions (23) we have 


L sln 2 7 „ 

D 2 - 



7 « 1 



p Q = 0.0027 slugs/ft 3 


p = 


23,500 


ft 


-1 


The choice Vi = 1 implies that the trajectories are those corresponding 
to entry from circular satellite orbit , with the altitude at entry being 
in the neighborhood of 80 miles (condition c). 


Nonlifting Vehicles 


Figure 1 shows the deceleration history in g r s for nonlifting 
vehicles having initial flight-path angles of 0°, 2°, and 4°. For all 
but a portion of the 4 case, the decelerations are seen to be below the 
level of 11.5 g considered to be an upper bound of human tolerance in 
reference 3* Hence, for the chosen set of initial conditions, the entry 
angle range 0 ° to 4° appears to cover the spectrum of usable trajectories 
for manned nonlifting vehicles. The case 7 j_ = 22° is included in fig- 
ure 1 to illustrate the magnitude of decelerations typical of trajectories 
having steeper entry angles, such as are experienced by long-range bal- 
listic missiles. This case will again be used in subsequent figures as 
illustrative of ballistic missile trajectories, in order to compare the 
oscillatory behavior of manned vehicles with that typical of ballistic 
missiles . 

Effect of K . - The effect of the dynamic stability parameter, K, on 
the oscillatory motion can be demonstrated by use of the convergence 
criterion, equation ( 38 ). Since the sign of the expression determines 
whether a vehicle T s oscillations are convergent or divergent, setting the 
expression equal to zero gives the set of circumstances signifying a 
changeover from convergent to divergent oscillations or vice-versa. Then, 
at each point on the known trajectory, one can solve for the value of K 
that makes the expression zero. The locus of such values of K plotted 
against altitude therefore forms a boundary which separates the range of 
altitudes over which oscillatory divergence is possible from the range 
over which it is not. Such boundary curves are shown in figure 2 for the 
three cases ( 7 ^ = 0 °, 2 °, 4°) considered suitable for manned vehicles and 
also for the ballistic missile case, 7 j_ = 22°. The drag-loading parameter 
W/CdA is 30 pounds per square foot . 2 

The significance of the figure is as follows: Vehicles whose values 

of K and altitude fall within the boundaries of their respective curves 

2 The assumption of constant aerodynamic coefficients is undoubtedly 
invalid at subsonic and transonic speeds. For this reason the curves are 
terminated below an altitude corresponding to sonic flight speed 
(V ~ 0.04). It would be possible to continue the curves if the varia- 
tions with Mach number of the aerodynamic coefficients were given; one 
would use equation ( 36 ) rather than ( 38 ) in this case. 
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will experience divergent oscillations over the altitude range within the 
boundaries. Thus, for example, the curve for the ballistic missile 
( 7 j_ = 22°) indicates that a missile whose value of K is -0.4 will 
experience divergent oscillations during its descent through the altitude 
range from 110,000 feet to 72,000 feet. On the other hand, for a missile 
with K - +0.4, the oscillations begin to diverge at 130,000 feet and 
continue to diverge thereafter. 

The main point that becomes evident from figure 2 is that the con- 
sequences of a deficiency of aerodynamic damping (K > 0) can be more 
serious for manned vehicles than they are for ballistic missiles. This 
follows from the observation that for the same value of positive K the 
amplitude of oscillation of the manned vehicle begins to diverge at a 
higher altitude, so that by the time the amplitude of the ballistic 
missile also begins to diverge the manned vehicle has already sustained 
a divergent oscillation over an altitude span of some 30,000 to 55*000 
feet. The reason for this becomes clear from inspection of the case 
K = 0, corresponding to zero aerodynamic damping. Referring to equa- 
tion (37) * we see that divergence begins in this case when q ! (s) = 0; 
that is, at the altitude where dynamic pressure is a maximum. Hence y 
any factor (smaller entrance angle being one) which tends to raise the 
altitude at which q is maximum will also tend to increase the serious- 
ness of the dynamic stability problem for cases where K > 0. 

Amplitude ratios . - The value of a plot such as figure 2 is that one 
can determine immediately whether or not the possibility exists of a 
dynamic stability problem for a vehicle with a given value of K. If a 
problem is indicated, however, its severity is still in doubt, since the 
amplitude ratio may be so small when divergence begins, or divergence may 
occur at so low an altitude, that the amplitude cannot grow to serious 
magnitudes in the time remaining before the vehicle lands . To investigate 
the severity of the problem, the actual amplitude history is required. 

Amplitude ratios have been evaluated from equation (27) for values 
of K of -2, 0, and +2, representative of vehicles with a large amount 
of aerodynamic damping, zero aerodynamic damping, and a large deficiency 
of aerodynamic damping, respectively. These are plotted as functions of 
altitude for a habitable vehicle (y^ = 0°) in figure 3(a) and for a bal- 
listic missile (y ± = 22°) in figure 3(b). The drag-loading parameter is 
again 30 lb/sq ft in both cases. In addition, the dimensionless velocity 
V is shown. 

The results show that with K = -2 the amplitude ratio diminishes 
very rapidly for both vehicles. For zero aerodynamic damping (K = 0), 
the motion is divergent below the altitude for maximum q, but the 
amplitude does not grow excessively for either the manned or ballistic 
vehicle before sonic velocity is reached. For K = +2, however, it is 
clear that a deficiency of aerodynamic damping of this magnitude leads 
to a divergent oscillation that may be of serious consequence. In this 
case the fact that divergence begins at a higher altitude for the manned 



1 6 


vehicle than it does for the ballistic missile makes the dynamic stability- 
problem for the former vehicle considerably- more severe. Thus, the ampli- 
tude ratio for the manned vehicle has already reached 1 at an altitude of 
115,000 feet (fig. 3(a)) whereas the amplitude ratio of the ballistic 
missile at this altitude is still only 0.13 (fig* 3 (b)). 

Effect of 7 j[ and W/CpA.- Having determined that the possibility of 
a serious dynamic stability problem exists only for positive values of K, 
we shall confine the remainder of the discussion to this case, letting 
K = +2. 

Figure 4 shows the effect on amplitude ratio of varying the initial 
flight-path angle for one value of W/CpA, whereas figure 5 shows 

the effect of varying w/CpA for one value of 71 . It will be seen that 
increasing either y± or w/CpA tends to diminish the severity of the 
dynamic stability problem for manned vehicles, although the former effect 
is small because of the small range of entry angles available for feasible 
manned flight trajectories. Both of these effects are simply manifestations 
of the single fact already pointed out, that factors tending to lower the 
altitude at which q is maximum also tend to diminish the severity of the 
dynamic stability problem. 

One final point should be made regarding the effects of 7 ^ and W/CpA 
on the oscillatory behavior of ballistic missiles and manned vehicles. 

To isolate the effect of entry angle in figure 4 the two vehicles were 
given the same value of W/C^A, namely, 30 lb/sq ft. We recognize that 
while this is a reasonable figure for the manned vehicle, the drag-loading 
parameter of a ballistic missile will generally be considerably larger. 
Hence, in general, the ballistic missile will benefit from larger values 
of both 7 i and W/CpA. Thus, in figure 5 even the curve for the manned 
vehicle with the highest value of W/CpA, 100 lb/sq ft, would fall above 
and to the right of a typical ballistic-missile curve. It is the combi- 
nation of beneficial factors, larger 71 and larger W/CpA, that leads 
one to the conclusion that for a ballistic missile, dynamic instability 
is a potential source of difficulty only in its terminal phase of flight. 

In contrast, it is seen that dynamic instability is a source of difficulty 
for the manned vehicle over a significant portion of its trajectory. We 
should note, however, that a compensating factor also exists for the 
manned vehicle that may ease the problem of controlling oscillatory diver- 
gence: The maximum value of dynamic pressure experienced by manned 

vehicles will be very much less than that of ballistic missiles. Hence, 
although control of the oscillatory divergence of the manned vehicle may 
be more urgently required, the frequencies at which the controlling device 
must operate, being proportional to will be very much smaller for the 

manned vehicle than they are for the ballistic missile. 



Lifting Vehicles 


To study the effect of lift, we have chosen a vehicle having a lift- 
to-drag ratio of 0.5 and a drag-loading parameter of 30 lb/sq ft. The 
range of initial flight-path angles is 0 ° to 4°, as for the nonlifting 
vehicles. The degree to which even this small amount of lift is effective 
in reducing decelerations may be illustrated by noting that the maximum 
deceleration experienced by the lifting vehicle in the worst case 
( 7 ^ = i|°) is only 3.6 g as compared to 13-3 g for the nonlifting vehicle 
with 7 j_ = 4°. Hence we are assured that at least from the standpoint of 
tolerable decelerations the vehicle under consideration is habitable. 

Effect of K . - A comparison of the convergence boundary curve for the 
lifting vehicle with the curve for the nonlifting vehicle is shown in fig- 
ure 6 for one entrance angle, y± = 0°. The figure shows that for vehicles 
deficient in aerodynamic damping (K > 0) the use of lift increases the 
severity of the dynamic stability problem in about the same proportion as 
it was increased by the use of small initial flight-path angles for non- 
lifting vehicles (fig. 2). Consideration of the K = 0 case reveals that 
this is once again attributable to the fact that the altitude at which q 
is a maximum has been raised still further. 

Amplitude ratios . - Amplitude ratios for the lifting vehicle with 
7 j_ = 0° and K = -2, 0, +2 are shown in figure 7 . These are to be compared 
with the results for the equivalent nonlifting vehicle given in fig- 
ure 3(a). It is noted that for K = 0 (and presumably for a range of 
small K) the amplitude history for the lifting vehicle is not signifi- 
cantly worse than that for the nonlifting vehicle. For K as large as 
+ 2 , however, it will be noted that the divergence is considerably more 
severe . 

Amplitude ratio and velocity histories for the lifting vehicle with 
initial flight-path angles of 2 ° and b° are shown in figures 8 (a) and 8 (b), 
respectively. For initial angles other than 0 , the trajectory of the 
lifting vehicle differs from that of the nonlifting vehicle in that the 
lifting vehicle makes small skips as the lift force momentarily overcomes 
gravity. The strange behavior of the amplitude ratio curves simply 
reflects the behavior of the static flight trajectory. Thus as the 
vehicle first gains and then loses altitude, the concurrent changes in 
dynamic pressure cause the motion first to diverge and then to converge. 
Finally, when the velocity has slowed sufficiently, the oscillatory 
amplitude behavior is almost identical to that for 7 ^ = 0 (fig. 7)* 


CONCLUDING REMARKS 


An analysis has been carried out of the oscillatory motions developed 
by vehicles as they traverse arbitrarily prescribed trajectories through 
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the atmosphere. Expressions were derived that describe the oscillatory 
motions and, in particular, the envelope of oscillations, as continuous 
functions of the properties of the trajectory. For the special case of 
constant aerodynamic coefficients the envelope of oscillations has the 
simple form 



where q and u are, respectively, the dynamic pressure and horizontal 
component of velocity, K is a measure of the aerodynamic damping proper- 
ties of the vehicle, and C is an arbitrary constant. 

Results of the analysis were used to study the oscillatory behavior 
of vehicles which have deceleration histories that remain within bounds 
of human tolerance. It was determined that for vehicles deficient in 
aerodynamic damping (K > 0), the factor governing the seriousness of the 
divergent oscillations which occur is the altitude at which dynamic pres- 
sure is a maximum. The decelerations of manned vehicles are maintained 
within tolerable limits by the use of small initial flight -path angles or 
by the use of lift. Since both tend to raise the altitude at which 
dynamic pressure is maximum, it was concluded that the dynamic stability 
problem of manned nonlifting vehicles deficient in aerodynamic damping 
may be more severe than that of comparable ballistic missiles, and that 
the use of lift further increases the severity of the problem. On the 
other hand, it is recognized that a compensating factor exists in that 
the maximum value of q experienced by manned vehicles will be very much 
sm all er than that of ballistic missiles. This factor may serve to ease 
the problem of controlling the oscillatory divergence of manned vehicles, 
since the frequencies at which a control must operate, being proportional 
to •R, will be very much lower for manned vehicles than they are for 
ballistic missiles. 


Ames Research Center 

National Aeronautics and Space Administration 
Moffett Field, Calif., Dec. 2, 1958 
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Figure 6.- Comparison of convergence -boundary curves for lifting and 
nonlifting vehicles with zero initial flight-path angles. 
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